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Observational evidence for AGN fueling. I. The merging of 

NGC 6104 with a companion. 

A. A. Smirnova, A. V. Moiseev, and V. L. Afanasiev 
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Abstract. We investigate in details the kinematics and morphology of the Seyfert galaxy NGC 6104 in order 
to identify the mechanism of gas transportation to the active galactic nucleus (AGN). Our observational data 
were obtained at the 6-m Special Astrophysical Observatory telescope with the MPFS integral-field spectrograph 
and the SCORPIO universal device in three modes: direct imaging, a scanning Fabry-Perot interferometer, and 
long-slit spectroscopy. Images from the HST archive were invoked to study the structure of the circumnuclear 
region. An analysis of deep images has revealed that NGC 6104 is in the phase of active merging with a companion 
galaxy. We have been able to study the detailed picture of ionized gas motions up to galactocentric distances of 
14 kpc and to construct the stellar velocity field for the inner region. The radial gas motions toward the AGN 
along the central bar play a significant role at galactocentric distances of 1-5 kpc. In addition, we have detected 
an outflow of ionized gas from the nucleus that presumably resulted from the intrusion of a radio jet into the 
ambient interstellar medium. Using diagnostic diagrams, we estimate the contributions from the AGN and star 
formation to the galactic gas ionization. We estimate the bar pattern speed by the Tremaine- Weinberg method 
and show that the inner ring observed in the galactic images has a resonant nature. Two possible ring formation 
scenarios (before and during the interaction with a companion) are discussed. 



1. Introduction. 

The problem of galactic activity still remains to be far 
from understood. Why are the nuclei active in some galax- 
ies and inactive in other ones? What is this, a short phase 
in the lifetime of any galaxy or a fundamental difference 
of the nuclei? There is no unequivocal answer to all these 
questions as yet. The mechanism of gas transportation to 
the region controlled by the gravity of the central object (a 
supermassive black hole, in the common opinion) remains 
a main problem. More specifically, much of the angular 
momentum must be drained from the gas that was ini- 
tially at galactocentric distances of several kpc to ensure 
its motion toward the nucleus (see the reviews of Combes 
2001; Wada 2004). Numerous attempts to associate the 
presence of an active (Seyfert) galactic nucleus (AGN) 
with properties of the host galaxy, such as the presence 
of a central bar or even a double bar (Laine et al. 2002; 
Knapen 2005) or an inner minispiral (Martini et al. 2003), 
have been made over the past decades. 

Many authors have tried to find a correlation between 
the presence of an AGN in a galaxy and its environ- 
ment: the presence of companions or traces of interac- 
tion (Dahari 1985; de Robertis et al. 1998; Schmitt 2001; 
Knapen 2005). However, statistically significant correla- 
tion of the properties of host galaxy with activity wasn't 
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found in any of the listed papers. Kauffmann et al. (2004) 
and Hall and Richards (2004) compiled various samples 
of galaxies with active and normal (quiescent) nuclei from 
the Sloan Digital Sky Survey and compared their morpho- 
logical properties. Again, no statistically significant differ- 
ences between galaxies with active and normal nuclei were 
found. 

In a recent paper devoted to this problem (Martini 
2004), it has been suggested that the mechanism of exci- 
tation and maintenance of activity could be complex and 
depend on several factors at once. The paper by Garcia- 
Burillo et al. (2005), in which the authors attempted to 
directly measure the molecular gas inflow rate into an 
AGN, is of considerable interest in this respect. A detailed 
analysis of the dynamics of inner regions in four Seyfert 
galaxies leads to a paradoxical conclusion - despite the 
presence of a bar, the gravitational potential distribution 
in these galaxies currently prevents the gas transportation 
into the innermost (less than 100-200 pc) region. Authors 
suggested several possible solutions of this paradox: al- 
lowance for viscosity and tidal friction that take away the 
angular momentum of gaseous clouds or a periodic action 
of some additional factors help the gas reach the nucleus 
(a mass redistribution in barred galaxies). In any case, 
it is clear that we can understand why the fuel reaches 
the region of gravitational influence of the central engine 
of active galaxies only through a detailed analysis of the 
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Table 1. Main characteristics of NGC 6104 



Type of activity (from NED) 


Syl.5 


Morphological type (from NED) 


S(R)Pec 


Morphological type (from LEDA) 


SBab 


Systemic velocity (from NED) 
Absolute magnitude 1 


8429 km/s 


-21.14 


Adopted distance 2 


113 Mpc 


Spatial scale 


0.55 kpc/" 



magnitude in the B-band was taken from 
Ho & Peng (2001) and given at distance 113 Mpc 
2 according to Mouri et al. (2002) 



kinematics and dynamics of inner (100-1000 pc) galactic 
regions. 

With this paper, we begin a series of papers devoted 
to observational evidences of gas inflow into the center, 
which must ensure the (possibly periodic) existence of an 
AGN. Integral-field spectroscopy allows us to trace the 
behavior of gas and stars both in the immediate vicinity of 
an AGN and far from the center. In addition, the method 
of ionization diagrams means to determination the type 
of an ionization source in various galactic regions. 

Here, using the galaxy NGC 6104 as an example, we 
wish to trace how the interaction between galaxies affects 
the activity of their nuclei. On the other hand, it is inter- 
esting to consider the relationship between nuclear activ- 
ity and peculiarities of the motion of gas and stars in the 
galactic disk. 

NGC 6104 (its basic parameters are given in Table 
1) is interesting in that no candidates for the role of a 
perturbing companion have ever been found (Rafanelli et 
al. 1995), although some of the authors pointed out its 
peculiar morphology (Jimcncz-Bcnito et al. 2000; Malkan 
et al. 1998; Tomita et al. 1999). Only our deep images 
have revealed outer asymmetric structures of low surface 
brightness that can be unambiguously interpreted as a re- 
sult of a recent interaction with a companion that was 
most likely followed by its disruption. Therefore, it is in- 
teresting to know whether the radial gas flows toward the 
galactic center that we detected are related to these signs 
of interaction, or not. 

2. Observations and data reduction. 

All of the observational data were obtained at the prime 
focus of the 6-m Special Astrophysical Observatory (SAO) 
telescope. The central region of the galaxy was observed 
with the multipupil fiber spectrograph (MPFS). In order 
to obtain information about the outer parts of the galaxy 
and its surroundings, SCORPIO focal reducer was used in 
the following observational modes: a scanning Fabry-Perot 
interferometer (FPI), long-slit spectroscopy, and direct 
broad-band imaging. The CCD detectors were TK1024 
(1024 x 1024) and EEV42-40 (2048 x 2048 pixels) in 2002 
and 2003-2005, respectively. A log of observations is given 
in Table 2. 



2.1. MPFS 

MPFS (Afanasiev et al. 2001) can simultaneously take 
spectra from 240 spatial elements arranged in the form of 
a 16 x 15 lens array with a scale of l"/element. The spec- 
tral range included both numerous emission lines of ion- 
ized gas ([OI]A6300, [OIII]A4959+5007, [NII]A6548+6583, 
[SII]A6717 + 6731, the Balmer H a and H^ lines) and ab- 
sorption features characteristic of the galactic old stel- 
lar population. We used a software package running in 
the IDL environment to reduce the observational data. 
The MPFS data reduction steps were briefly described by 
Moiseev et al. (2004). The data reduction result is a data 
cube in which a spectrum corresponds to each image pixel. 
NGC 6104 was observed twice in the same spectral range. 
After the primary data reduction, the two cubes were com- 
bined. The resulting field of view of the cube was 22 x 20" 
, while the spatial resolution after the smoothing proce- 
dures corresponded to a seeing of <~ 2". The spectra of 
spectrophotometric standard stars were used in order to 
transform the galaxys spectra to energy units. 

The stellar velocity fields were constructed using a 
cross-correlation technique modified to work with MPFS 
data (Moiseev 2001). Because of the relatively low spec- 
tral resolution, we could not study the distribution of the 
stellar velocity dispersion within the field of view and can 
only conclude that it does not exceed 200 km s -1 . 

By fitting the profiles of the main emission lines with 
Gaussians, we constructed their surface brightness distri- 
butions and radial velocity fields. The maps for the most 
intense lines are shown in Figs. 1 and 2. The influence of 
stellar absorption features, primarily in the spectral re- 
gion with hydrogen lines, constitutes a main problem in 
analyzing the spectral lines of ionized gas. We used the 
following simple technique that allows the contribution 
from the stellar continuum to be taken into account, at 
least to a first approximation, without using complex and 
not always unambiguous methods like synthetic spectra. 
The absorption features under the H Q and Hp lines were 
fitted with Gaussians, the positions of whose centers were 
determined from the stellar velocity field (strictly speak- 
ing, from its fit by the model of circular rotation) , FWHM 
was fixed (because of the low spectral resolution, we may 
ignore the stellar velocity dispersion variations over the 
field), and the amplitude was proportional to the contin- 
uum level 1 (i.e., corresponded to a constant equivalent 
width). The free Gaussian parameters (the absorption line 
FWHM and equivalent widths) were chosen by fitting the 
spectra at the edge of the field of view, where the contribu- 
tion from the emission lines decreases and the absorption 
features under H a and H3 can be clearly detected. 



1 To take into account the contribution from the AGN- 
related non-thermal component, we subtracted the image of 
a point source from the continuum maps in the innermost re- 
gions. The resulting maps corresponding to the stellar compo- 
nent alone were used for calculation of the equivalent widths. 
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Table 2. Log of observations NGC 6104 and its environment 



Date 


Instrument 


Exposition 


spectral 


spectral 


seeing 






(sec) 


range 


resolution 




2002 May 15 


MPFS 


3 x 900 


4600 - 7000A 


8A 


1.5" 


2002 Aug. 13 


MPFS 


2 x 1200 


4600 - 7000A 


8A 


1.5" 


2003 Dec. 24 


SCORPIO (FPI) 


32 x 90 


H a 


2.5A 


1.6" 


2004 Jan. 31 


SCORPIO 


3 x 300 


B 




1.6" 




(direct images) 


3 x 180 


V 




1.6" 






8 x 120 


R 




1.5" 


2004 June 21 


SCORPIO 


4 x 300 


4000 - 9000A 


10A 


1.1" 




(long slit) 










2005 May 17 


SCORPIO 


2 x 900 


4020 - 5600A 


5A 


1.6" 



(long slit) 



2.2. SCORPIO focal reducer. 

SCORPIO universal instrument (Afanasiev and Moiseev 
2005) allows various types of spectroscopic and photomet- 
ric observations to be performed with ~ 6' field of view. 
Below, we describe in details each of the modes used. 

2.2.1. Fabry-Perot interferometer. 

The observations with the scanning FPI were performed 
in the H a emission line; the spectral interval between the 
neighboring interference orders was ~ 28 A. The required 
spectral range was separated out by a narrow-band filter 
with FWHM of ~ 15 A. During the observations, we se- 
quentially took 32 interferograms of the object by chang- 
ing the FPI plate gap. The data reduction technique was 
described by Moiseev (2002a). The data reduction result 
is a data cube in which a 32-channel spectrum with a sam- 
pling step of 0.9A corresponds to each pixel with the scale 
lpix=0'/71. The radial velocity field and the monochro- 
matic H Q image were constructed by fitting the emission 
line profile with a Gaussian. Final spatial resolution af- 
ter the smoothing procedures corresponds to a seeing of 
2'.'1. The Haflux map was calibrated in absolute energy 
units (erg s" 1 cm~ 2 acrsec~ 2 ) by comparison with a sim- 
ilar MPFS map for the galactic center. 

2.2.2. Direct imaging. 

We obtained images of NGC 6104 in the Johnson-Causins 
B, V, R c bands with a scale of O'.'35/pixel. The data re- 
duction was performed in IDL and included a number of 
standard procedures: bias subtraction, flat fielding, and 
cosmic ray particle hits removal. The magnitude calibra- 
tion was performed using secondary standards. For this 
purpose, we obtained short-exposure images of both the 
galaxy itself and fields with photometric standards us- 
ing the Zeiss-1000 SAO telescope. Having calibrated these 
images using a standard technique, we were able to per- 
form the photometric calibration of our deep images. The 
achieved surface brightness levels (for a signal-to-noise ra- 
tio of 3) were 26.5, 25.4, and 25.5 m /D" for the B, V, and 
R c images, respectively. 



2.2.3. Long-slit spectroscopy. 

To study the galaxys environment, we took long-slit spec- 
tra of four faint galaxies, candidate dwarf companions 
of NGC 6104 (every time the slit crossed two neighbor- 
ing objects) and a spectrum of the asymmetric feature in 
the image of the elliptical galaxy CGCG 196-022 close to 
it. The data were reduced using a software package run- 
ning in IDL. The reduction steps were briefly described 
by Afanasiev and Moiseev (2005). 

3. Surface photometry for NGC 6104. 

Figure 3 shows (a) the galactic image from the Digital 
Palomar Observatory Sky Survey, (b) our deep R c -band 
image, and (c) the WFPC2/HST image. The HST ob- 
servational data have already been briefly described by 
Malkan et al. (1998), who pointed out an asymmetric "per- 
turbed" shape of the tightly wound spiral arms emerging 
from the bar. However, an examination of all the available 
(including H Q ) images leads us to conclude that a pseudo- 
ring is most likely present here, while the asymmetric spi- 
rals are traceable in the outer regions at r > 20". 

Figure 4 shows the radial variations in the shape of 
elliptical isophotes with radius (r). The cllipticity peak 
at r > 10" corresponds to the end of the circumnuclear 
bar elongated along the position angle PA = 88°. As r 
increases, the isophotal ellipticity decreases sharply, while 
PA remains almost constant up to r « 15". This behav- 
ior is attributable to the above mentioned ring with a 
slight inner ellipticity that is noticeable in all images and 
that has bluer colors than the rest of the galactic disk. 
This feature caused by the presence of a younger stellar 
population, the emission from HII regions ionized by O-B 
stars is the most expressive indicator of ongoing violent 
star formation in the ring. Thus, according to the FPI 
map, 70% of all galactic H Q luminosity is concentrated 
in the ring (and only 5% in the AGN). The total ob- 
served H Q luminosity of NGC 6104 is 2.5 x 10 41 er ffS - 1 
or 1.3 x 10 42 erg s" 1 after its correction for the internal 
extinction in the galaxy 2 . According to the formula from 

2 Here, we assumed that the internal extinction in all star- 
forming regions is the same, A(H a ) » 1.8 m (the reddening in 
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(a.) 



8200 8300 8400 8500 8600 




Fig. 1. Radial velocity fields of NGC 6104 (in kms^ 1 ): (a) 
in H Q emission line, constructed from FPI data; MPFS ob- 
servational data in the (b) H Q , (c) [Nil] , (d) [OIII] emis- 
sion lines and (e) the stellar velocities. The isovelocities 
are everywhere drawn at intervals of 30 kms -1 , the cen- 
tral contour corresponds to the assumed systemic velocity 
of 8418kms- 1 . 



Kennicutt (1998), a relatively high star formation rate of 
SFR w 10 Mayr- 1 , with 7 ' M Q yr^ 1 for the ring alone, 
corresponds to this luminosity. 

H a that we obtained from the H Q -to-H^ flux ratio using MPFS 
spectra for the HII region northeast of the nucleus). 



The shape of the ring deviates from a circle even in the 
galactic plane. This is clearly shown by Fig. 3c, which was 
constructed according to the disk orientation obtained by 
analyzing the gas velocity field (see below). We see that 
the ring remains elliptical (with an axial ratio of ~ 0.85); 
the boundary ring radii were assumed to be 11" and 18". 
The southern half of this structure forms an almost regular 
semiellipse, while its northern half is more irregular (per- 
turbed) in shape. The almost rectilinear fragment of the 
spiral arm between the bar and the inner edge of the ring 
looks particularly strange. Also, the ring center formally 
determined from isophotes is displaced by several arcsec- 
onds to the southeast of the nucleus. The displacement of 
the center of the star formation ring is also immediately 
apparent in the H a image (Fig. 2a). Below, we will return 
to the discussion of peculiarities in the shape of this struc- 
ture. 

The outer shells firstly found by us in the galaxy are 
its most intriguing morphological features. These shells 
are completely invisible in the images from the Digital 
Palomar Observatory Sky Survey (Fig. 3a) and are barely 
detected in the images from the Sloan Digital Sky Survey 
(SDSS). However, they are clearly distinguishable in our 
images in all three color bands at a surface brightness 
level of 24.5 — 25.5 m /n". These shells are most conspicu- 
ous southwest of the galactic center at distances of more 
than 50" (28 kpc) from the nucleus, i.e., in the opposite 
direction and considerably farther than the northeastern 
feature in the surface brightness profile (Fig. 3b). Many 
authors who classified the galaxy as a peculiar one have al- 
ready called attention to this feature (see the Introduction 
for references). Interestingly, despite the presence of this 
asymmetric feature in the images, the mean brightness 
profile (averaged in elliptical apertures elongated in the 
agreement with the galactic disk orientation) is satisfac- 
torily described by a classical exponential law. When the 
radial brightness profiles were examined, it was found that 
the spatial resolution of the images from the 6-m tele- 
scope (~l'/5) was too low to reliably separate the contri- 
butions from the compact bulge and the starlike nucleus. 
Therefore, we used an image through the F606W filter 
from the HST archive for a detailed analysis. 

We decomposed the surface brightness profile into 
standard components: an exponential disk (the radial 
scale length rdisk =9" 4 coincided with the scale deter- 
mined from the deep R c -band image), a de Vaucouleurs 
bulge (with an effective radius of r e =1"7), and a star- 
like nucleus. The derived component luminosity ratio 
Lbuige/ Ldisk = 0.14 corresponds to the Sc morphological 
type (according to Simien and de Vaucouleurs, 1986). 

After the deprojection of the images for NGC 6104 
to the face-on position and the subtraction of the two 
dimensional dusk+bulge model constructed from the sur- 
face brightness profile decomposition, it became clear that 
the outer shells, which are clearly seen in the images in all 
three color bands, are a single structure that bends around 
the disk of NGC 6104 (Fig. 3d). This is most closely resem- 
bles a companion galaxy that has been disrupted almost 
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Fig. 2. Monochromatic images of NGC 6104: (a) H Q (FPI) 
map. MPFS maps: in the (b) H Q , (c) [Nil], (d) [OIII] lines 
and (e) in the continuum near A6600A. 

totally by tidal forces. Similar structures that are formed 
when galaxies are disrupted by tidal forces have been de- 
scribed in many papers, as well observational ones as de- 
voted to numerical simulations (e.g., Jogee et al. 1999). In 
the case under consideration, we most likely observe the 
merging of a companion with NGC 6104. 

According to our estimates, the total luminosity of the 
outer low-surface-brightncss filament accounts for ~ 2% 
of the galactic total R-band luminosity (minus the lu- 
minosity from the starlike AGN). The fractional lumi- 



nosity increases to 3-7% if the bright fragments located 
to the northwest and the northeast of the nucleus at 
r = 20" — 30" are assumed to also belong to this tidal 
structure (Fig. 3d). Thus, there can only be merging ex- 
actly with a dwarf companion (with Mr > — 17 m ) whose 
mass accounts for no more than 2-7% of the mass of 
NGC 6104. 

4. Spectrophotometry and ionization diagrams. 

The galactic nucleus, bar and a bright field star are dis- 
tinguished in the MPFS continuum images. The bar and 
the star formation ring are clearly seen on the maps in the 
most intense emission lines (except [OIII] A5007). In the 
[OIII] image, only the bright nucleus is seen and there are 
no features far from the center. 

To determine the ionization sources of the central 
galactic regions, we constructed ionization diagnostic di- 
agrams from MPFS data. Using the ratios of lines with 
different excitation mechanisms, we can determine the 
regions where thermal (young hot stars), non-thermal 
(AGN), or shock ionization dominates. The boundaries 
that separate the regions with different ionization sources 
(AGN, HI or LINER) in Figs. 5a and 5b were taken from 
Veilleux (2002). We distinguished three regions in the 
MPFS field of view: a nucleus, a bright region at the end 
of the bar in the Ha, [Nil] and [SII] lines and points be- 
longs to the bar at distances of 3" — 7" from the nucleus. 
In both our diagrams, the points belonging to the nucleus 
lie in the region of ionization by nonthermal radiation. 
However, according to the [OIII] /Hp flux ratio criterion, 
about one third of all points from the nuclear region lie 
near the AGN-LINER boundary, i.e., the gas can also be 
partially ionized here by shocks. In the next section, we 
will return to this issue in the discussion about nuclear 
jet. 

All points from the region of the northeastern knot fall 
into the photoionization region; on our H a map (Fig. 2a), 
one of the brightest HII regions at the end of the bar corre- 
sponds to these points. In both diagrams, the points from 
the bar region are located in the transition region between 
AGN ionization and shock ionization; several points in the 
[SII]/H a diagram fall into the region of photoionization by 
young hot stars. Ionization by shocks, which can be asso- 
ciated both with the edges of the bar (where dust lanes 
are noticeable in the HST image) and with the nuclear 
jet (which, however, may be not so large, see below), is 
apparently also added to the two main competing mecha- 
nisms (photoionization/non-thermal source) at the points 
belongs to the bar region and the AGN. 

5. Analysis of the kinematics of ionized gas and 
stars. 

We determined the orientation parameters of the galactic 
gaseous disk using the model of pure circular rotation of 
a thin disk in which the inclination to the line of sight, io, 
and line of nodes position angle, PA , did not change with 
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Fig. 3. (a) POSS2 image of the NGC 6104 field, the numbers mark four candidates for dwarf companions, (b) R c image 
of NGC 6104 in the magnitude scale, the outer isophote corresponds to a surface brightness of 25 m /n", the isophotal 
step is l m . (c) HST image with the F606W filter deprojected to the face-on position. The dotted line corresponds to 
the line of nodes, the ellipses mark the adopted maximum and minimum radii of the ring with allowance made for its 
cllipticity. (d) Residual brightnesses of the R c -band image for NGC 6104 after deprojection to the face-on position 
and the subtraction of the disk model. 



radius. Here, we analyzed the large-scale H Q velocity field 
(excluded the points at r < 9", where significant peculiar 
motions take place). The following values were obtained by 
X 2 minimization of the observed velocities deviations from 
the model approximation: PAq = 56 ±3° and io = 39 ±5°; 



these values are in a good agreement with the data of the 
orientation of the outer disk isophotes (Fig. 4). 

Then, we analyzed the noncircular motions by the 
tilted-ring method. The velocity fields were broken down 
into elliptical l"-widc rings oriented according to the val- 
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Fig. 4. Results of the isophotal analysis of the images for 
NGC 6104 in three color bands: radial variations in the 
position angle of the major axis of the isophotes (a) and 
their ellipticity (b). The dotted line marks the location of 
the line of nodes. 



ues of PA and io obtained above. In each ring we de- 
termined the optimal position angle of the kinematic axis 
PAkin, the mean rotation velocity V rot and the systemic 
velocity V sys in the approximation of circular rotation. 
See, Moiseev et al. (2004) for details of the method and 
references to original papers. The center of the rings was 
fixed at the galactic nucleus. The symmetry center of the 
velocity field is slightly (by l'/5) displaced relative to the 
isophotal center of the continuum image, we attribute this 
displacement to the effect of the noncircular motions that 
we detected. The constructed dependences PAfc,„(r) and 
V rot (r) are shown in Figs. 6a and 6b. The velocity fields 
of ionized gas and stars are in a good agreement with one 
another, demonstrating a similar turn PAkin relative to 
the line of nodes (PAq) in the bar region (Fig. 6a). The 
relatively large scatter of MPFS data points is a result of 
a lower MPFS measurement accuracy than the accuracy 
of kinematic estimates in H Q with FPI, which has a higher 
spectral resolution and a smaller pixel size. However, there 
are also differences in the kinematics of gas and stars. 

Firstly, Fig. 6b shows that the rotation amplitude for 
the stellar component at r = 3" — 8" is appreciably (by 
~ 50 kms -1 ) smaller than for the ionized gas in H Q . This 
mismatch results from a considerably higher stellar veloc- 
ity dispersion than the gas one (the asymmetric drift effect 
widely known in stellar dynamics). More unusually, about 
half of all velocity measurements in the [Nil] doublet are 
also indicative of a systematically lower rotation velocity 
than that in H Q . In our opinion, this is attributable to 
a gas deceleration at the shock front that emerge at the 
bar edges where the fraction of the gas emission in forbid- 
den lines, including that in [Nil], increases. This effect was 



considered by Moiseev (2002b) for several barred galaxies. 
Recall that, according to the MPFS diagnostic diagrams, 
the emission from some regions inside the bar just points 
to the shock ionization mechanism. 

Secondly, an excess of "blue" ionized gas radial veloc- 
ities, which is particularly clearly seen in Fig. la, is ob- 
served in the nucleus. As a result, the formally measured 
V sys at r < 3" in the ionized gas velocity fields are ap- 
preciably lower than the mean systemic velocity at large 
radii (V sys = 8418 kms -1 ), while in the stellar component 
V S y S {r) ~ const. We interpret this excess of negative ve- 
locities as a signature of radial gas motions toward the ob- 
server; the outflow source itself is not resolved (i.e., its an- 
gular size is less than 1.5—2"). We observed similar signs of 
nuclear outflow when analyzed integral-field spectroscopy 
for the galaxies Mrk 315 (Ciroi et al. 2005) and NGC2273 
(Moiseev et al. 2004). Note that in all these objects, we 
observe not the high- velocity jet itself emerging from the 
AGN, but the result of the interaction between the radio 
jet and interstellar clouds, the formation of an expanding 
cocoon of hot gas. One of the most characteristic exam- 
ples was considered by Capetti et al. (1999), who observed 
Mrk 3. The size of such a structure (radio jet+cocoon) in 
Seyfert galaxies is only several hundred parsecs and it is 
not surprising that our spatial resolution is not enough to 
study it in NGC 6104 in detail. In NGC 6104, our measure- 
ments give an estimation of the gas outflow velocity along 
the line of sight - 30 - 40 kms" 1 in the H Q , [Nil], [SII] 
and [OIII] lines. In the latter case, the formal application 
of the tilted-ring model yields PAkin ~ 80°, which differs 
fundamentally from the measurements in the remaining 
lines. Here, most of the [OIII] emission apparently origi- 
nates from the circumnuclear jet and it is inappropriate 
to use the model of circular motion. 

It is possible that the points corresponding to the nu- 
cleus in the ionization diagrams (Figs. 5a and 5b) lie near 
the AGN/LINER boundary owing to the interaction of 
the jet with the interstellar medium. It should be noted 
that we give only a lower limit for the outflow velocity, be- 
cause the spectral resolution is too low for the component 
associated with the gas surrounding the jet to be directly 
identified; we detect only the shift of the line profile as a 
whole. Finally, note that we have failed to find any refer- 
ences in the literature to radio observations of NGC 6104 
with an angular resolution high enough to detect this ra- 
dio jet. 

PAkin within 10" of the center turns in the direction 
opposite to the line of nodes compared to PA of the bar 
isophotes (Fig. 6a). Such behavior of the position angle is 
in agreement with both available observational data on the 
kinematics of barred galaxies and with numerous model 
calculations of radial motions toward the center in the re- 
gion of a triaxial bar potential (for references to original 
papers, see Moiseev et al. 2004). This turn of PAkin is bet- 
ter seen in the FPI data. Because of the above-mentioned 
radial outflow from the nucleus, we should be very guard- 
edly in the estimates of PAkin at r = 1" — 3" based on 
such a simple rotation model, since the velocities of two 
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Fig. 5. (a, b) Diagnostic diagrams for NGC 6104. The filled circles, diamonds, and crosses, correspond to the nucleus, 
the star-forming region in H Q , and the bar region, respectively, (c) The mask for separating the above-mentioned 
regions with H a isophotes superimposed on it. 



spatially separated components, the gaseous disk and the 
outflow from the AGN, are projected onto the line of sight. 

Significant radial (inflow) motions take place only in- 
side the bar, since the kinematic axis coincides with the 
line of nodes of the disk even at r > 10" (i.e., in the ring 
region). We tried to estimate the amplitude of the radial 
motions using a simple model that assumed each point at 
a given radius r to follow circular rotation (V v ) and to 
be involved in radial motions (Vr), so the observed radial 
velocity is the sum of the components of both vectors: 

V bs{r) = V ays + V R (r) sin ip sin i + V v (r)cosip sin i. (1) 

Here, the azimuthal angle in the galactic plane <p is mea- 
sured from the line of nodes. We obtained Vr = 50 ± 
lOkms" 1 at r = 2" - 8" and V R w at r > 10", in 
agreement with the previous conclusion. 

Since few direct observational measurements for in- 
flow motions in bars are founded in the literature (they 
are considerably fewer than the bar pattern speed esti- 
mates), each such measurement is of considerable interest 
in discussing the role of bars in AGNs fueling and secular 
evolution of galaxies. Our value (~ 50kms _1 ) is close to 
that obtained by Jogcc ct al. (1999) for the nuclear bar of 
NGC 2782 (15-45 kms -1 ), but is considerably lower than 
the velocity of inflow motions in the much higher-contrast 
bar such as in NGC1530 (85 - 200 kms" 1 ; Regan et al. 
1997). 

6. Origin of the ring structure. 

As we have already noted above, a ring ~ 6 kpc in radius 
that differs from the remaining galactic regions by a high 
star formation rate is clearly seen in all optical images of 
NGC 6104. The following scenarios of the formation of 
this structure are possible. 

The first, most obvious cause for a barred galaxy is the 
formation of a ring at one of the dynamical resonances. 
Comparison of the results of numerical simulations with 



real ring galaxies (for references, see Buta and Combes 
1996) shows that the rings are formed near the inner and 
outer Lindblad resonances (the nuclear and outer rings, re- 
spectively) and at the ultraharmonic (1:4) resonance near 
the end of the bar (the "inner ring" ) . The size and loca- 
tion of the ring in NGC 6104 exactly corresponds to the 
latter case. 

The second cause is considerably rarer, but possible 
in an interacting system is a collisional ring, i.e., an ex- 
panding ring density wave generated by the passage of a 
massive companion through the galactic disk. However, 
we have found no ring expansion in NGC 6104; signifi- 
cant non-circular motions (different from an expansion) 
are observed only in the inner region associated with the 
bar. Similarly, for obvious reasons that follow from the 
galactic morphology and kinematics, a polar ring is not 
possible either. 

A further argument for the resonant nature of the ring 
is its elongated shape with an axial ratio of 0.85 in the 
galactic plane. The ring is elongated in the same direction 
as the bar, which must be in the case for rings at res- 
onances within the corotation radius, in agreement with 
observations (Buta 1986) and model calculations (Schwarz 
1984). Therefore, it would be interesting to independently 
determine the positions of the main resonances in the 
galaxy by measuring the bar pattern speed Q p . 

Numerous methods for estimating £l p are described in 
the literature (see Hernandez et al. 2005), but the over- 
whelming majority of them are indirect ones based on 
additional assumptions about the bar formation mech- 
anisms or on numerical simulations. However, in recent 
years, a model-independent method based on considerably 
more general postulates formulated twenty years ago by 
Trcmaine and Weinberg (1984) has come into wide use. 
Let a rotating (gaseous or stellar) medium distributed in 
a flat disk, whose motion is characterized by a single pat- 
tern speed fl p and for which the continuity equation holds, 
in any case, on time scales of several disk rotations. The 
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observed parameters at each point in the disk are its sur- 
face brightness and line-of-sight velocity. Then, 

(V) = (X)n p smi, (2) 

where the x -axis is parallel to the galaxy line of nodes, 
(V) is the weighted mean observed radial velocity in a nar- 
row strip along this axis, and (X) is the weighted mean 
coordinate along this strip; the surface brightness is used 
as a weight factor. An additional condition is a linear re- 
lation between the surface brightness and surface density 
of the medium under study; the latter must become zero 
at the disk edges. The Tremaine- Weinberg (TW) method 
implies an infinitely thin, flat disk and ignores any verti- 
cal motions. Note that, despite the traces of interaction 
in the galactic outer regions, the inner gaseous disk of 
NGC 6104 may be considered flat, at least where there 
are points in the FPI velocity field. Indeed, as follows 
from Fig. 6a, PAkin is almost equal to the mean value 
at r = 10" — 25"; moreover, there are no significant vari- 
ations in the parameters of elliptical isophotes at larger 
distances, r = 35" - 55" (19-30 kpc) (Fig. 4). Therefore, 
the warp of the disk is possible only in the outermost re- 
gion. 

The TW method was used for measuring fi p in the 
stellar components of SBO galaxies (Aguerri et al. 2003; 
Debattista 2003) and in the molecular disks of late-type 
galaxies (Zimmer et al. 2004; Rand and Wallin 2004). 
Here, we do not consider the possibility that the valid- 
ity conditions of the TW method are satisfied in concrete 
objects as well as the possible error sources, referring the 
reader to the above-mentioned papers. 

Hernandez et al. (2005) applied the TW method to 
the FPI ionized gas velocity field. However, in our point 
of view, the reasoning of these authors regarding the appli- 
cability of the TW method to the ionized gas distribution 
is not convincing, primarily because there is no simple 
relation between the gas density and its surface bright- 
ness in H Q . In the case of NGC 6104, however, it has 
already been pointed out that the gas and stellar velocity 
fields arc similar, the difference lies mainly only in the ro- 
tation amplitude, while the behavior of PAkin along the 
radius is the same. Therefore, to a first approximation, the 
H Q velocity field may be assumed to be equivalent to the 
stellar velocity field and, if we use the continuum bright- 
ness as a weight factor, then we can try to apply the TW 
method to our FPI observations. Of course, such a some- 
what liberal treatment of the TW method requires more 
careful reasoning, but preliminary results show that the 
values of Q p obtained by the method described above for 
several galaxies observed with the 6-m telescope are in sat- 
isfactory agreement with other methods for estimating the 
positions of resonances (Afanasiev and Poltorak, private 
communication) . Thorough substantiation and testing the 
applicability of the TW method to these data require an 
individual study that is outside the frame of this paper. 

We took several strips parallel to the line of nodes 
through the velocity field and the galaxy continuum im- 
age and calculated (X) and (V) in each of them. The strip 



width and step was 2"1, i.e., 3 pixels. Figure 6c shows a 
clear correlation between {X) and (V), as follows from Eq. 
(2). Linear fitting of this dependence yields the bar pat- 
tern speed Ct p = (19±4)kms _1 kpc _1 (to within the sign 
defined by the chosen direction of the x axis). Here, the 
uncertainty in choosing the position angle was assumed to 
make the main contribution to the error in the O p estima- 
tion (Debattista 2003). We varied PA within PA ± 5°, 
which, in general, slightly exceeds its formal measurement 
error, and determined how much Q p changed in this case. 

Figure 6d shows the radial dependence of the angular 
velocity f2 and the ± n/2 and f2 — k/4 curves, where 
k is the epicyclic frequency. In calculating Q and k, we 
used a seventh-degree polynomial fit to the galactic rota- 
tion curve to avoid any drops in the derivative dQ/dr. An 
examination of the constructed resonance diagram shows 
that the adopted bar pattern speed intersects the O — k/4 
curve at a distance of 17±3". That is the 1:4 resonance co- 
incides with the ring, within the measurement error limits, 
as follows from the theory. Moreover, at this bar pattern 
speed, the fl p — const line does not intersect the il — k/2 
curve. In other words, the bar hasn't inner Lindblad res- 
onance (ILR) where the radial gas motions toward the 
center cease (Combes and Gerin 1985), so gaseous clouds 
can reach the AGN, providing fueling for it. 

It should be emphasized that the main conclusion that 
follows from Fig. 6d, the absence of ILR, can also be 
reached without using the TW method. It will sufficient 
to assume that the ring is at the 1:4 resonance and this 
is the most likely assumption, as we noted at the begin- 
ning of this section. It then follows from the diagram of 
characteristic frequencies that the bar pattern speed is 
O p = 17 — 27kms~ 1 kpc -1 , which matches the results of 
the TW method, within the error limits. 

7. Examination of the galactic environment. 

Based on images from the 6-m telescope, we have detected 
asymmetric filaments of low surface brightness in the outer 
regions of NGC 6104 located at apparent galactocentric 
distances of 25-45 kpc. In contrast to spiral density waves, 
these structures are clearly non-stationary and must dis- 
sipate in a time of the order of several periods of disk 
rotations at a given radius. Under the assumption of a 
flat rotation curve, we find the characteristic lifetime of 
the filaments is 0.6 — 1.0 Gyr, which imposes a constraint 
on the maximum age of the event that generated them. As 
we have already noted above, these structures resemble in 
shape the result of the disruption of a companion galaxy 
by tidal forces. Moreover, they are very similar to the fil- 
aments near the disk galaxy NGC 2782. In the latter ob- 
ject, radio data on the HI distribution convincingly prove 
that a merger with a companion with a mass of ~ 1/4 of 
the galactic mass took place (Jogee et al. 1999). Besides, 
after the deprojection of the NGC 6104 image to the face- 
on position and the subtraction of a two-dimensional disk 
model, it became clear that the outer shells, which are 
clearly seen in the images in all three color bands, are 
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Fig. 6. Analysis of the inner kinematics of NGC 6104. (a) Radial variations in the position angle of the kinematic axis. 
The dotted line indicates the orientation of the line of nodes and the solid line indicates the location of the central 
bar. (b) The observed rotation curve of the galaxy; the solid line indicates a polynomial fit. (c) Determination of the 
bar pattern speed by the Tremaine- Weinberg method, (d) The diagram of characteristic frequencies. The gray color 
marks the range of possible bar pattern speeds (f2 p ) the dotted lines indicate the maximal and minimal ring radii. 



a single structure warps around the disk of NGC 6104, 
confirming once again the hypothesis of a disrupted com- 
panion. In general, however, an alternative explanation is 
also possible - these structures are tidal in nature and re- 
sulted from a perturbation of the galactic disk during a 
close passage of a companion without its disruption. The 
companion must then be still close enough. Let us consider 
this case in more detail. 

NGC 6104 belongs to the cluster Zwicky 1615.8+3505, 
also known as N45-389 in the X-ray field. According to 
Tomita et al. (1999), this is a dynamically young poor 
cluster 4 Mpc in diameter. Although the cluster is fairly 
poor, the contribution from the possible interactions of 
NGC 6104 with its other members cannot be disregarded. 

Four extended objects that can play the role of dwarf 
companions are seen in the POSS2 image in the immediate 
vicinity of NGC 6104 (within three galactic optical diam- 
eters). These are marked by numbers in Fig. 3a. Since 
the distances to them are unknown, we obtained spec- 
tra for all four objects. Analysis of the spectra containing 
both absorption and emission features shows that these 
are relatively distant galaxies (with redshifts of 0.096, 
0.220, 0.316, and 0.097, respectively) and their proxim- 
ity to NGC 6104 is a projection effect. 



Another candidate companion is the elliptical galaxy 
CGCG 196-022 (CGCG 1615.0+3550), at 4' from 
NGC 6104, and, according to Tomita et al. (1999), its 
luminosity is a factor of three lower. The radial veloc- 
ity difference between NGC 6104 and CGCG 196-022 is 
205 kms -1 (from NED). Since the velocity dispersion of 
the cluster galaxies is high (584 kms -1 in Tomita et al. 
1999), we estimate a minimal time of possible encounter 
between these two galaxies, 0.2 Gyr ago, which is con- 
siderably shorter than the characteristic lifetime of the 
filaments. Thus, if this galaxy had a specific trajectory, 
it could encounter NGC 6104 and produce the mentioned 
above tidal structures; but the galaxy CGCG 196-022 it- 
self must also be perturbed appreciably, since the masses 
of the two galaxies are similar. However, analysis of our 
deep images obtained with the 6-m telescope for CGCG 
196-022 has revealed no perturbations of its isophotes, ex- 
cept an elongated feature on the periphery of the galaxy 
southwest of the nucleus, which is also clearly seen in Fig. 
3a. We obtain a spectrum of this feature. In the range 
4020 - 5600A (see Table 2), we detected only one, but 
fairly seen intense emission feature at a wavelength of 
4362A. For moderate redshifts, such situation is possible 
only if we observe the [OII]A3727 line; thus, the veloc- 
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ity of this structure is ~ 51200 kms -1 and it is a distant 
background edge-on galaxy. 

Thus, although the only one candidate galaxy that 
could affect appreciably the gravitational potential of the 
disk in NGC 6104 is CGCG 196-022, there is no conclu- 
sive evidence for their recent interaction. On the other 
hand, since CGCG 196-022 is an elliptical galaxy, the 
tidal distortions in it may be considerably smaller than 
those in NGC 6104, a late-type disk galaxy. However, the 
appearance of the external structures more closely corre- 
sponds to the disruption of a dwarf companion (a single 
strongly warped structure; see above for the analogy with 
NGC 2782) than to the tidal tail from the interaction with 
a galaxy of comparable mass. In the latter case, one might 
expect the presence of a bridge between the galaxies. This 
is further confirmation that the outer filaments in the im- 
age of NGC 6104 are the result of the disruption of a 
companion rather than the tidal perturbation in the disk. 

8. Conclusions. 

We investigated in detail the active galaxy NGC 6104 
and its environment by means of the integral-field spec- 
troscopy. We constructed the ionized gas and stellar ve- 
locity fields and studied the contributions from various 
mechanisms to the gas ionization in the inner (r < 5 kpc) 
region using diagnostic diagrams. Based on surface pho- 
tometry, we refined the morphological type of the galaxy, 
Sc. During our work, we were able to analyze the pecu- 
liarities of the structure of NGC 6104 on various spatial 
scales. 

1. In the innermost region of the galaxy, we suspected 
the existence of a jet from the AGN with a size of no 
more than several hundred parsecs. To be more pre- 
cise, we found effects that could be explained in terms 
of the interaction between a radio jet and the ambi- 
ent interstellar medium: an excess of "blue" velocities 
in all lines of ionized gas (up to 40kms _1 along the 
line of sight), a peculiar appearance of the [OIII] ve- 
locity field, and the possible contribution from shocks 
to the ionization of the ambient gas that follows from 
diagnostic diagrams. To our regret, there are no radio 
observations of NGC 6104 with an angular resolution 
high enough to detect this radio jet. 

2. At larger galactocentric distances (1-5 kpc), we identi- 
fied radial motions of ionized gas along the bar toward 
the nucleus with velocities of ~ 50kms _1 . Diagnostic 
diagrams and the observed deceleration of rotation in 
the [Nil] line relative to H Q point to an appreciable 
contribution from the shock fronts at the bar edges to 
the ambient gas ionization. We measured the bar pat- 
tern speed by the Tremaine- Weinberg method. Since 
the inner Lindblad resonance was shown to be absent 
in the bar, the radial gas motions can reach the AGN, 
fueling its central engine. 

3. We considered the circumnuclear ring and showed that 
70% of all the current star formation in the galaxy 



(SFR w 7M Q yr~ 1 ) proceeds in it. The ring is at ul- 
traharmonic (1:4) resonance of the bar. We noted a 
distorted shape of the ring (first of all, a displacement 
of the center relative to the nucleus) ; it could be that 
this asymmetry was produced by a recent interaction 
with a relatively massive companion. 
4. Finally, in the outermost regions of the galaxy at galac- 
tocentric distances of 25-45 kpc, we detected extended 
asymmetric filaments of low surface brightness that are 
clearly seen in all three bands, BVR C . Similar shells, 
the so-called ripples, can be formed in galaxies during 
a close passage of a companion galaxy or during ac- 
cretion of dwarf galaxies (Schweizer and Seitzer 1988; 
Kormcndy 1989). Analyzing the images from the 6-m 
telescope, we found in the residual images (after the 
subtraction of a disk model) that these are not two 
different shells, but a single structure that resembles 
most closely a companion disrupted by tidal forces. 

So, is a relation between the large-scale outer struc- 
tures (the remnants of the merged satellite) and such in- 
ner features of the galaxy as the radial motions toward 
the center and existence of an AGN? Here, two alterna- 
tives are possible: either the ring and the bar existed in 
the galactic disk before the interaction with a compan- 
ion or they emerged already during this interaction. In 
any case, the interaction leads to a mass redistribution in 
the galaxy, primarily to a higher mass concentration in 
the central (one kiloparsec size) region (Combes 2001). If 
the bar and the ring already existed in the galaxy, then 
such a redistribution could increase the rate of gas inflow 
through the bar into the AGN. A number of peculiarities 
in the shape of the ring noted above (asymmetry, displace- 
ment of the center) can then be explained by distortions 
of the gravitational potential via the interaction. In this 
case, the merger with a companion must have taken place 
much later than ~ 1.0 Gyr determined in the previous 
section, since the dynamical times at r = 5 kpc are an 
order of magnitude shorter than those at r = 50 kpc and 
the ring had not yet been perturbed and relaxed. 

The more courageous assumption that both the bar 
and the ring are the result of interaction with a compan- 
ion is supported by the following circumstantial evidence. 
We see from Fig. 3c that the concentration of dust lanes 
in the HST image is highest in the northwestern half of 
the galaxy. If the dust lies in the disk plane, then this half 
must be closer to the observer. Since it follows from the 
velocity field that the northeastern edge of the galaxy re- 
cedes from us, then we find that the galactic disk rotates 
counterclockwise. However, the spirals in the disk are then 
the leading ones; this is a very rare situation that is possi- 
ble only in interacting systems in which the spiral density 
wave is generated by a companion (Pasha 1985). We then 
find that the entire spirals-ring-bar structure resulted from 
the gravitational tidal effect of a companion (the situation 
in which the bar and the spirals rotate in opposite direc- 
tions seems a short lived one and breaks up after several 
bar rotations). This scheme is very attractive, but note 
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that the dust distribution in the circumnuclear region is 
still very complex and further arguments for such a spatial 
orientation of the galactic disk are required. 

Note once again a number of analogies between 
NGC 6104 and NGC 2782, which was studied in detail 
by Jogee et al. (1999). Both objects are disk spiral barred 
galaxies with their outer isophotes perturbed by a recent 
interaction that ended with the merging with a less mas- 
sive companion. In both cases, outer filaments with a sim- 
ilar structure (the remnants of the disrupted companion) 
are observed. In NGC 2782 this interaction led to the con- 
centration of gas in the inner region, the formation of 
a bar from molecular gas, and a powerful starburst. In 
NGC 6104, we also see gas flows toward the center where 
the AGN is located. Although it is not yet clear why simi- 
lar causes led to the formation of different types of nuclear 
activity (a Seyfert nucleus and a circumnuclear starburst), 
a relationship with a recent merging event is traceable in 
both objects. 

According to the studies listed in the Introduction, 
which are based on samples of both Seyfert and normal 
galaxies, no statistically significant difference has been 
found between these two classes in terms of the presence 
of companions or traces of interaction. In most studies, 
the authors restricted themselves to searching for traces 
of large-scale interactions with companions of comparable 
mass (major- merging). However, a more detailed study of 
the surroundings of galaxies can revise these conclusions; 
the apparently single (at first glance) galaxy can be in 
the process of merging with a low-mass companion, as we 
found in our case. A similar situation is also observed, for 
example, in the active galaxy Mrk 315. Previously, it was 
considered in all papers as a single galaxy, but deep im- 
ages from the 6-m telescope (Ciroi et al. 2005) show that 
this galaxy simultaneously interacts with two dwarf com- 
panions, one of which is in the process of tidal disruption. 
Finally, we must emphasize that the AGN activity in any 
galaxy is not necessarily caused by the presence of a bar 
or traces of a recent interaction. We try to draw attention 
to the fact that only a careful analysis of each object us- 
ing currently available observational facilities allows the 
factors that caused its current activity to be established. 
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